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^observations extended, it is easy to calculate approximately the 
j|correction required by Leverrier’s value of the parallax. As¬ 
suming the average tabular parallax of Mars during the Ascension 
^observations to have been 21", and taking Mr. Grill’s constant 
w ^ ’376 as the correction to be deducted from this value, we have 
for the correction to the Sun’s tabular parallax of 8"‘95, 

*•376 x 8-95 _ ff , l6 
21’00 

Thus, we get 8 // *95 — 7/ *i6 — 8 // *79, agreeing extremely well with 
the general result of Mr. Grill’s observations. 

I have for some time thought that an extension of the 
method employed by Leverrier in the case of the Sun to deter¬ 
mine the parallax of some of the nearer planets, would probably 
be attended with success, and the result brought out incidentally 
by Mr. Grill seems to support this view. One of the main diffi¬ 
culties would be the accurate determination of the Moon’s mass. 
Probably, sufficiently numerous and exact observations exist to 
determine the lunar parallax to within tooVoo^ of its actual 
value, i.e. to about "*03, while the length of the seconds pen¬ 
dulum may be considered as already known with an equal 
exactitude; these two data should therefore give us the value of 
the Moon’s mass to within the yoVo^ 1 part of the truth. As 
respects other methods of determining the solar parallax, it 
appears to me that something more might probably be accom¬ 
plished by a careful analysis of existing observations of Jupiter's 
satellites, treating the diameters of the planet and of the satellites 
as unknown quantities to be determined from the eclipse obser¬ 
vations, and taking into account the possible effect of the refrac¬ 
tion of the Jovian atmosphere (which may be considerable) on 
the times of disappearance and reappearance. 

Stroud Green, N., 

1881, July 13. - 


On the Determination of the Moon's Libration. 

By A. Marth, Esq. 

In the method which has been adopted in several investiga¬ 
tions of the Moon’s libration, the apparent position of a spot 
referred to the centre of the Moon’s disk and to the circle of 
declination is deduced from the measured distances and position- 
angles of a number of points of the Moon’s limb in reference to 
the spot; from the apparent coordinates thus found the seleno- 
centric right ascension and declination of the spot are deduced, 
and these are converted into ecliptical longitude and latitude. 
Comparison is then made with the selenocentric longitude and 
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IJatitude computed from approximate values of the spot’s seleno- 
■jgraphical coordinates and of the inclination of the lunar equator 
;3jo the ecliptic, and the equations are found which exhibit the 
Inflect of corrections of these values and of the terms due to 
“physical libration. 

It appears preferable to follow another course; to proceed 
from assumed selenographical coordinates of the spot to its 
apparent coordinates referred to the centre of the Moon’s disk 
and the direction of the Moon’s axis, to determine the correc¬ 
tions of these coordinates which best satisfy the measurements, 
and to form the equations between the corrections and the effect 
of the Moon’s physical libration, &c., upon these coordinates. 

Let I denote the assumed inclination of the plane of the 
lunar equator to the ecliptic; & + 180 0 the longitude of its 
ascending node, which coincides with the descending node of the 
orbit; J its inclination to the plane of the Earth’s equator; N" 
the right ascension of the ascending node. If, further, l Q is the 
mean longitude of the Moon in its orbit and M+180 0 the arc of 
the lunar equator between equator and ecliptic, let l 0 — 4-M-fhf 
be denoted by a c . The departure of the first or zero meridian, 
from which the selenographical longitudes are reckoned, from 
the ascending node & +180° of the lunar equator on the ecliptic, 
is, in consequence of the equality of the periods of rotation and 
mean revolution, = Z 0 —■■ Q>, and therefore its departure from the 
ascending node of the lunar equator on the Earth’s equator 


= l Q — 8 + M + i8o° = a 0 —N + i8o°, 


so that the angle u 0 is in reference to the equator analogous to 
l 0 in reference to the ecliptic. 

Ephemerides of the position of the lunar equator referred to 
that of the Earth have been given in the Berliner Jahrbuch since 
1839, and in the Nautical Almanac since 1867, the notation and 
formulae adopted being those employed by Encke in his paper 
u Ueber die selenocentrischen Constanten bei den Sternbe- 
deckungen und die Berechnung der Libration des Mondes, 
nebst Tafeln,” in the Berliner Jahrbuch for 1843. Instead of 
the arc of the lunar equator between the Earth’s equator and 
the ecliptic, therein given, it will be found far more convenient 
to have the value of a OJ which is in Encke’s notation 

*= Iq — S3 +18o° + A -f~ £3 t 

For the present purpose the Ephemerides will have to be 
computed rather more accurately, l Q and £3 being duly referred 
to the true equinox. Since, moreover, other formulae will be 
found preferable to Encke’s, there would be no advantage in 
retaining his designations i, A and £V, which are replaced by 
the letters J, M+180° and N in the present paper. 

The values of a 0 , hf and J may be deduced from those of l Q > 
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! J? , I, and of tlie apparent obliquity e of the ecliptic by means of 
lilhe formula 

Is! tan 7 } = -tan I cos 


■co 1 
loo I 

I'—11 


tan N == — tan I sin Q> 
tan J = 
sin J (a 0 -l 0 ) = 


COS 7 ] 

‘ sin (e + rff 

tan (e + n) sec N, 

sin A I sin 1 e . _ 

- 2 - 2— . sin Q . 


COS -g J 


The next proceeding is the computation of the position-angle 
P of the Moon’s axis, and of the selenographieal longitude A and 
latitude B of the centre of the Moon’s disk as seen by the ob¬ 
server. Let a f and c f be the apparent right ascension and decli¬ 
nation of the Moon’s centre, which are derived from the geocentric 
right ascension and declination a and c>, if 0 is the sidereal time, by 

r' cos S' sin (V — a) = — p' sin (0 — a), 

r f cos S' cos (ct! — a) =r cos S — p' cos (0 — a), 
r 1 sin 5' = r sin S — p", 

or by equivalent formulae, in which r(=— - 1 - — ) and r' are the 

linear distances of the Moon’s centre from the Earth’s centre 
and from the observer, and p' and p u the perpendicular distances 
of the observer from the Earth’s axis and from the plane of the 
equator, the unit of these distances being the equatoreal semi¬ 
diameter of the Earth. The selenocentrie spherical triangle 
between the pole of the lunar equator, the pole of the celestial 
equator and the observer furnishes the equations : 


cos 13 sin P = — sin J cos (a'—N), 
cos B cos P= + sin J sin (a' —N) sin 5' +cos J cos S', 

sin B = -f sin J sin (a' —1ST) cos S' —cos J sin S', 

cos B sin (A + a Q — N) = cos J sin (a' —N) cos 5' + sin J sin S', 

cos B cos (A + a 0 — N) = cos (a' — N) cos S'; 

from which latter two may be deduced 


cos B sin (A + a Q — a') = [sin S' — tan \ J sin (a' — N) cos 5'] sin J cos (a' —N) 
= (sin S' — sic B) tan J J cos (a' —N). 

Hence are derived the formulae for the computation of P, B 
and A:— 

tan £= tan J sin (a'— N), 

tan P = — tan J cos (a' —N) —-- , 
v cos (£—S ) 

tan B = tan (f — 5') . cos P, 

• /. , sin B —sin 6' . -r> 

sm A + « c -a = --■ — t - • Bin P, 

2 COS 2 i J 

sin -A (B — Sf) cos \ (B + S') . -n 

or =-2A- > —LA- L . sm P, 

cos 2 ^ J 

or — . sin P , if tan f r = tan A J sin (a'—N). 

cos C 2 ^ 1 
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Then follows the computation of the apparent position of the 
■|jpot. If X and j 3 are the selenographical longitude and latitude 
the spot, and p its position-angle at the centre of the disk, as 
jljeen by the observer, we have selenocentrically, if the angular 
“distance between spot and observer is called S, 


sin S sin (p — P) = cos j 8 sin (A — A), 

sin S cos (p — P)= —cos jS cos (A —A) sin B- 8 cos B, 

COS S = COS j3 COS (A —A ; cos B - - in 8 sin B. 


If the linear distance of the spot from the Moon’s centre of 
gravity is k r and from the observer r" , the unit of distances being, 
as before, the semi-diameter of the Earth’s equator, the angular 
distance s of the spot from the centre of the disk, as seen by the 
observer, and r" are found by 

r" sin 8 = Jc' sin S, 
r" cos s=r'—Jc' cos S; 

or, indirectly, by 

Jc 

sin s = — sin (S + s), 
r' 

. sin S 

r ' = r' . -•. 

sin (fS + s) 

The coordinates 

sc" = co" sin s sin (p — P), (w" =*= —I— = 2 o 6264 " , 8 ), 

arc 1 

y ,f = sin s cos (p — P), 


or the apparent rectangular coordinates of the spot in reference 
to the centre of the disk and to the computed direction of the 
Moon’s axis seem to me the most suitable for the determination 
of the libration. The corrections of x” and y" which best 
satisfy the observations, compared with the variations of x" and 
y" which represent the effect of the corrections of the assumed 
values of I, X, / 3 , and of the different terms of the physical libra¬ 
tion, furnish then the equations of condition for the determination 
of the latter. 

In the differentiation of the equation for x n and y ,f 

„ w'W • c • , n\ to "Jc l n . 

w rf =——. sm S sin (p—P) = —— . cos j8 sm (A—A), 

f^n yfr 

03” ¥ a''k f 

p" = ——. sin S cos (p — P) = _ 7 _[ — cos £ cos (A—A) sin B + sin 0 cos B], 
it must be borne in mind that P as a function of the assumed 
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.^elements is itself variable, so that if the corrected values x" -f lx n 
!|hnd y" + ty" are to refer to the same direction P as x n and y n 
gthemselves, the variation of P must be duly taken into account, 
jrapifferentiating accordingly, we get 


8x"=~X 


— ^ cos i 8 cos (A - A) 8 (A — A) 
— sin /3 sin (A —A) 8 $ 

+ sin S cos (jo— P) 8 PJ. 


, „ 8k f 

* 



+y" - 


cos $ sin B sin (A —A) 8 (A —A) 

+ [cos 0 cos B + sin j 8 sin B cos (A—A)] 8/3 

- [cos /3 cos B cos (A - A) + sin £ sin B] SB 

— sin S sin (p — P) SP X 

W J 

V * 


Prom the expressions of A, B, P are derived 

SA = - 8a Q - [tan B sin (a 0 -JST + A)~ tan 1 J] sin J.SN 
— tan B cos (a 0 — K + A) SJ, 

SB= -cos (a 0 -N + A) sin J. 8 N +sin (a c -N + A) SJ, 
cos B.5P= - sin (a 0 -N + A) sin J. 8 N-cos (a 0 -N + A) 8 J. 

But, if 

M = a 0 -N-(l 0 ~a), 

then 


8a 0 = 8l 0 - (tan |J cos M + tan | I) sin I.8Q 

— tan ^ J sin M.8I, 

sin J.8N = _ cos M sin 1.8 Q, — sin M.8I, 

= + sin M sin 1.8 S 3 — cos M.8I. 


Hence 


8 A = - « 0 + [tan B sin (l Q -Q + A) + tan J1] sin 1.8a 
+ tan B cos (l Q -Q +a) SI, 

5B = cos (I q —Q+ A) s in 1.8& - sin (l Q — Q + A) 81, 
cos B.5P = sin {l Q - & + A) sin 1.8 Q + cos (l Q - £ + A) 81, 
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;|he substitution of wbicb gives 

_ eo S (A-A) (S; o + SX) 

loo i / L 

i^—11 


— tan ft sin (A —A) 5j3 

+ [tan ft sin (l 0 — 9> + A) + tan % I cos (A — A)] sin 1.5 & 
+ tan ft cos ( l Q — Q> + A) 51 




425 


*y" ^ ^ cos B | - tanB sin (A —A) {§l 0 + 5A) 

+ [1 + tan ft tan E cos (A —A)] 5/8 
— [cos (l 0 — Q + A) + tan ft tan E cos ( l 0 — Q> + A) 

— tan J I tan B sin (A —A)] sin 1.5 & 




+ [sin (l Q — S3 + A) + tan ft tan E sin 


ft,-a+ a)] si} 


or, putting 


then 


1 + 


tan ft tan E sin (A — X)-v sin (A, — A) 
tan ft tan B cos (A — X) = v cos (A! —A), 


5y" = ^ c ^f — cos B <^ — tan B sin (A —A) (dl Q 4 - 5a) 

+ v cos (Aj — A) 5j8 

— \y cos (l Q — 9> + Aj) — tan | I tan B sin (A —A)] sin 1,5^ 


„ M f 

+ y-T- 


+ v sin {l Q — Q> + A,) 51 


The physical libration produces oscillations in the values of l Q9 
Q, and I, the theoretical investigation of which is to be found 
in the Mec. Gel. L. v., and in the papers of Poisson (Conn, des 
Temps, for 1821 and 1822), Bessel (Astr. JSfachr. Nos. 376 and 377, 
or Abhandlungen, iii. pp. 317-328), Wichmann (Astr. Nachr. Nos. 
619, 621, 628, 630, 631), and Hartwig ( Beitrag zur Bestimmung 
der physischen Libration des Mondes , 1880). If i denotes the 
inclination of the Moon’s orbit to the ecliptic, e its eccentricity, 
m the angle in the orbit from the node to the perigee, g the 
Moon’s mean anomaly and 2 (H sin A) the sum of the several 
inequalities of the motion in longitude, which produce sensible 
terms ; if, further, m!, n', 7 v r denote the daily motions of the mean 
longitude, of the node and of the perigee, and h' the daily motion 
of A, all of which quantities are known : the oscillations of Z 0 , & 
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Le to the physical libration, may be represented by the 
equations :— 

«o = 2(-H. 6 • sin A) + 0 sin (A 0 + v Q t), 

sin 15^ = sin I U 0 

— Kq cos o) sin (l 0 — f<!) 

+ nqf sin co cos (l Q — U ) 

+ b Q sin [B 0 +/ 4 ^— ft> — ft)] 

+ e Q sin (C 0 + fi't) cos (l Q — ft) 

+ cos (C 0 + fi't) sin (l Q — ft), 

V f 

51= — Kq cos co cos (l Q — ft) 

— Kqf sin co sin (l Q — ft) 

~b Q cos [3 0 + nt-(l 0 -ti)] 

— c Q sin (O 0 + p't) sin (l Q — ft) 

On 

+ — 7 = 2 =- cos (C Q + j u!t) cos (l 0 — Q), 

V f ■ 


^Ind I, dn 
following 


in which/, a oi A 0 , b Q , B 0 , c 0 , 0 o are the seven constants of the 
problem which are to be determined, and in which 6 , Kq, Kq f , v Q , 
jj, fj,' have the following values :— 


e' = 

e = 


— 2 7 vm 1 

h'h f ' I 7 l 

0 ' 


for each of the 
effective inequalities, 


1 - 0 '’ 

f, , 0-0391 n 

w' — n' 1 i + I )' 


. — 2n’ei 

K Q ~—r 


K —TV 


, - a/ _ 2 n ' m '. J — • a/i -/ 


/jl — rn! — - 


w'l 


i + I 


/ 4 ?/I — 

— • V/. 


3 i + I 


In order to compute the coefficients of the unknown quanti¬ 
ties correctly, the values of I and/employed must be sufficiently 
near the true ones, or the computation will have to be repeated 
with corrected values. Instead of / itself, it will be preferable 
to make the correction cf of the best available value of /, the 
unknown quantity, or rather 0/ multiplied by a properly chosen 
factor k", so that in the equation for cl 0 the term 


2(1 


H 0(i + 0) 

(1 +vy 


sm 


*)«"*/ 
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i to be added to the term. 


4 2 7 


which is to be considered as known. 

The substitution of the values of cL, c£l, and cl, which 
represent the effect of the physical librat:- u. in the previously 
given equations for c;e' and cy", furnishes then the equations: 

*. „ kJl' COS & 

Sx = — — X 

y cos (A — A.) — tan 0 tan I sin (l a — S, + A)^ ^( — — s * n 

+ q tan | I sin g cos (A—A) 

+ q tan j8 cos (cy + A) 

+ (q~qf) tan j3 sin cy sin A 

kJc' cos j 8 f , f A . v 8 a 
- ~ { + cos (A —A) . “ 

+ tan 0 sin (A—A). — 

K 

— tan j8 cos (Z Q — & + A) . -51 

ft 

+ £ [cos (A - A) - tan j8 tan I sin (/ Q - 8 + A)] ~ • ^ "^y r sin ^ 

— — tan /3 sin u> sin A I a"<5/ 

k" J 

+ ^cos (A — A) — tan /3 tan I sin (l 0 — 9> + A)^ , sin (A 0 + v Q t ) 

+ tan i3 . cos (B 0 + fit + A) 


K 

2 cos 


— tan j8 . -—— . 

Vf cos 


. c s cos (G 0 + fi f t-A l )\ } 
20S Aj K J 


if tan A, 


= V f 


tan A 


x' f ioo 


ioo 


V 


and, putting as before 

rsinf/q — A) = + tan 0 tan B sin (A — A), 
v cos (A,—A) = i + tan 0 tan B cos (A — A), 
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S' s ft kJc' cos j8 t-, 

S| fy" =-. COS B . If X 


IOOI 

"{( 


rp ! * 

tan I cos (l 0 — Q + A)+ - aD ^ sm 




+ q tan J-1 sin g . — — — sin (A — A) 
+ q sin (to + Aj) 

+ (q — qf) sin oo cos A, J> 


kJc' cos $ 
r n 


. cos B 


■’{ 


, tan B . %n 5a 
+ -sm (A — A)— 

V K 


+ COS (A,—A) 


5B 


— sin (^ 0 —^ + Aj) .— 

K 

+ [(tan I co S (^-fi + A) + ^ s in(A- A) ) sin j) 

■+■ -h- sin oo cos A,]*"?/- 

+ (tan I cos (l Q —£2 + A) + sin (A—A)\ . a °- sin (A 0 + vj.) 

\ VjK 

+ ^ sin (B 0 + fit + Aj) 

K 


+ 


cos Aj 

COS A, 


. ^sin (Co + ^-A^V 

K J 


+ 


y" ioo Ek' 


if tan A ;/ = —tan \ 

V f 


ioo 




The factor k, which is introduced on account of obvious con¬ 
venience, is arbitrary, and may be put =240 or, perhaps better, 
=300. The most convenient value of a", the factor of 3 /, will 
be that of the largest of the coefficients 

H 0(1 4 -0) 

* ’ (1 +fef 

or that dependent on the Moon’s annual equation. 

I have written out the equations for <ta?" and %y n in full, 
retaining all the terms in symbolical form, so as to leave the 
question of the best numerical values which ought to be 
employed an open one. The substitution of numerical terms 
and the introduction of suitable abbreviations will, of course, 
greatly simplify the equations, especially in case the selected spot 
is near the lunar equator, a case the advantages of which are 
apparent at a glance. 

The measures from which the observed values of lx n and 
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fy" are to be deduced furnish a series of angular distances 
Ijj . . . s n of the spot from points of the Moon’s limb in various 
osition-angles yq . . . the corresponding times being t x . . . t n , 
'or these times the values of 


a', 5', /, R', P, B, A, y" 


(R' being the Moon’s apparent semi-diameter, sin R' 



be found by interpolation from tbe values of these quantities 
computed at intervals of, say, 20 or 30 minutes for the extent of 
each set of measurements. From the distances and position- 
angles, duly cleared from instrumental errors and from the 
effects of refraction, the corresponding semi-diameters R x . . . R ri 
are to be derived by means of the equations 


s f sin (p,—P,)+ #/'** R sin Q /? 
cos fe-F/) +y/' = h cos Q„ 


s n sin (p n - P w ) + x n " = E n sin Q , 
S n COS (p n - P ri ) + Vn = R n COS Q n , 


in which Q| . . . Q n are the position-angles of the observed 
points of the limb at the centre of the disk, reckoned from the 
direction of the Moon’s axis. 

The comparison of these Rj . . . R % derived from observa¬ 
tion, with the interpolated values R'j . . . R' w supplies the 
equations of condition 

K/ — K - sin Q ; fa" + cos Q, 5y ”- ^ SR 0 , 

Eq 


R' n — R w = sin Q, n + cos Q w Sy" — 5_8R 0 , 

R 0 

from which the corrections dx", ly" and c)R 0 for each set of mea¬ 
sures are to be deduced. The correction c)R 0 is that of some 
conveniently chosen mean value R 0 , and, in case the tabular 
parallaxes can be depended on, is common to all the sets of 
measures made with the same instrument under similar cir¬ 
cumstances. 

The substitution of the corrections lx n and ly n derived from 
observation in the previously given equations leads to the equa¬ 
tions of condition for the determination of eft, M, If and the 
others terms of the physical libration. 

Tn order to be able to trace the influence of the irregularities 
of the Moon’s limb, it will be well to compute the seleno- 
graphical positions of the several points of the limb, which 
have been measured in position-angle and distance in reference 
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\j$o the spot. Their selenographical longitudes \ n and latitudes 
||J n are derived from the corresponding values of Q nf A n9 B n bj 
;5neans of the equations 

loo i 

cos J8 n sin (A n — A n ) = — sin Q n , 


cos j8 w cos (A„ — A n ) = - cos Q n sin B n + tan R n cos B n , 
sin p n = + cos Qjn cos ~B n + tan R n sin B n , 


which in practice, where great accuracy is not required, may be 
sufficiently simplified. 

Although I have left the question of the best numerical 
values which ought to be substituted in the formulae an open 
one, some remarks which bear upon it and touch a few of 
the practical difficulties of the investigation will not be out of 
place. 

An uncertainty in the determination of I is occasioned by the 
difficulty of disentangling £1 and b Q , in case the observations do 
not extend over a sufficiently long period. The angle B 0 -f pt is 

n 'I , 

= B 0 + ^- T — T Jf, 

2 + 1 

= -B 0 -(l 0 -U ) 0 + l 0 -Q + ^L.t, 

2 + -L 

if ( 7 0 — &3) 0 is the value of l Q — £3 for the time t= o, from which 
t is reckoned ; or, putting 


= and 
n'i 

B / + jut = B^ + l Q — Q + fi t t. 


Hence in the equation for ly" 


b Q sin (B 0 + (d + \) = (b 0 sin B,) cos ( l Q - Q, + A, + /*/) 
+ ( b Q cos B ; ) sin (l 0 - Q + A, + fi t). 


But the coefficient of —H is sin (Z c —S^+Xj), so that the 
two angles differ only by u { t. As the period of this difference is 
about twenty-four years, it is obvious that observations extend¬ 
ing only over a year or two do not allow the satisfactory deter¬ 
mination of both unknown quantities BI and b 0 cos B /5 and that 
one of the two must remain practically indeterminate. 

A troublesome difficulty arises from the necessity of repeating 
a portion of the computations, in case the adopted value / is not 
sufficiently near the true one. Wichmann’s first value/--0*2603 
was derived from a preliminary solution of his longitude- 
equations, and it served, together with Nicollet’s value of 
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1 

||=i° 28' 47 ", in the further computations. Wichmann’s final 
||leducfcions gave 

| I=I °3 2 ' 9"; /=g? = o- 4 i9, 

or, neglecting a Q and b 0 on account of the uncertainty of their 
determination, 

T o ' r / * S^S"-6 ^ _ 

1=1° 32 24 ; t ~--\ ; =0-42^. 

J ^ 1^4 

It is to be regretted that Hartwig, instead of adopting one of 
these latter values, has retained in his computations the pre¬ 
liminary inaccurate approximations. His own final result 
{Monthly Notices for May, p. 376), 

I=i° 3 6 ' 39 "; 7 = 0507 , 

gives for v Q the value or the period 888’6 days. Wich¬ 

mann’s (last) result gives 1520"'9, or period 852-1 days, while 
the preliminary value employed is 1728"'9. or the period 749’6 
days. How far the determination of a Q is affected by the error 
of the assumed value of v 0 is at present unknown; hut it need 
scarcely be more than hinted that the error affects most seriously 
the connection between the results of the observations of 1845 
and of 1878. Bessel points out in Article 6 of his paper that, 
in case a sensible value of a c should be successfully ascertained, 
the determination of the duration of its period would afford the 
surest means for the exact determination of f. According to the 
data on page 376 the value of a Q is found from 


Wichmann’s Observations = 183-2 ± 80 6 
Hartwig’s Observations = 62-3 ±41-9. 

It is to be hoped that before long the publication of the 
investigation of the Moon’s libration which has been under¬ 
taken at the Oxford University Observatory will facilitate some 
decision of what are the best available values of I and/. 

In a second paper I shall treat the problem of determining the 
physical libration by observations of the changes which it pro¬ 
duces in the relative position of spots ; I shall, further, show the 
effect of the libration upon the observations of a spot in right 
ascension and declination, and also consider some other questions. 
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